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Results from the next to leading order (NLO) dimuon analysis from the NuTeV experiment 
at FNAL are presented. Charged current interactions in neutrino-nucleon scattering with 
two oppositely charged muons in the final state allow direct study of charm production and 
measurement of the strange sea. NuTeV's sign selected beam gives it the ability to extract 
the strange and antistrange seas independently, for which an asymmetry has been predicted 
— , in some theoretical models, and which is currently of intense interest in interpreting neutrino 

£T) • electroweak results. The dimuon results presented here represent the first analysis of dimuon 

events performed utilizing DISCO, a fully NLO cross section code differential in all variables 
, required to model detector acceptance. 
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XI 1 Introduction 
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Events containing two oppositely signed muons (dimuons) from muon neutrino deep inelastic 
scattering experiments provide a unique window into the strange quark content of the nucleoffl 
These events occur in charmed particle production from a charged current (CC) interaction 
^ . with a strange (or down) quark. Approximately 10% of the time the charmed hadrons decay 

semi-muonically, yielding a final state with two oppositely charged muons, one from the weak 
vertex, one from the charm decay. Charm production from down quarks is Cabibbo suppressed, 
making dimuons most sensitive to the strange sea. Dimuons are clearly distinguishable in a large 
neutrino detector such as that employed by NuTeV^. 

NuTeV ran during FNAL's 96-97 fixed target run and recorded 5102 dimuon events from CC 
z/^, and 1458 from CC scattering in its iron target. The detector was calibrated throughout 
the run by muon, electron, and hadron beams so its response is well understood. NuTeV's 
beamline was constructed to be able to select or beams with high purity. This a priori 
knowledge of whether events were the result of a neutrino or antineutrino interaction, unique to 
NuTeV, allows one to measure the strange sea independently from the antistrange sea. 

Of particular topical interest is whether the strange sea is different from the antistrange 
sea. It has been proposed, for example, that the proton wavefunction could have a virtual 
K + K pair component which would lead to an asymm etry b etween strange and antistrange 
sea!3. More recently this possibility has been entertained ^^ as an explanation for the almost 
3d difference between NuTeV's sin 2 6w result and global fits. To eliminate this discrepancy a 
strange-antistrange asymmetry, as defined by the momentum weighted integral: S~ = J x[s(x) — 
s(x)]dx, would need to be positive and as large as +0.00',^. NuTeV's dimuon data sample offers 
the opportunity to address this question directly. 



2 Fitting the Dimuon Cross Section 



Nu TeV has made its dimuon data available in the form of a forward dimuon cross section 
tablPl. The cross section of dimuons with charm decay muons of energy greater than 5 GeV 
was extracted in bins of Bjorken x, inelasticity y, and neutrino energy E v , and corrected for 
detector smearing effects. 

A fit to the dimuon cross section table is made up of the following components: 
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Where is the measured dimuon cross section table mentioned above, to which the left 

side of the expression above is compared. da ^y m is the cross section model, dependent on 
the strange sea and charm mass, which would be varied in a fit. EMC is an x dependent 
correction for the EMC effect and nuclear shadowing, parameterized from a fi t to data from 
charged lepton scattering experiments in heavy nuclear targets and deuteriunU^l B n is the 
semileptonic charm branching ratio (0.093 ± 0.008, from a re- analysis of FNAL E-53£E1), and 
A is a kinematic acceptance correction accounting for the 5 GeV cut on the energy of the muon 
from the semileptonic charm decay. A depends on E u , x, y as well as fragmentation and, at 
NLO, charm mass. A Monte Carlo simulation using the DISC(PI 

NLO cross section code, 
differential in all variables required to model this acceptance, was used to calculate A. 

Fits were performed using acceptance corrections calculated with a constant Collins-Spillel^ 
epsilon of 0.75, since the charm mass and epsilon are correlated at NLO. It was determined that 
the strange asymmetry was insensitive to the choice of epsilon, and a value where NLO Monte 
Carlo matched to data well was chosen, and varied over a conservative range (0.25) for error 
estimates. The charm mass is then allowed to vary in the fit. 



3 Strange Asymmetry Results 

Several extractions of the strange and antistrange seas from the dimuon cross section tables have 
been performed at different orders in QCD (LO and NLO) as well as with different methods 
of parameterization. The first class of fits treat the strange and antistrange seas completely 
independently, and use a more or less traditional parameterization: 
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QCD evolution is approximated here by modifying the shape of already evolved u and d 
distributions through k, a, a. Results from a LO fit based on the CTEQ5lP PDF set are 
shown in the left half of figure ^ This and all subsequent plots are shown for a typical Q 2 of 10 
GeV 2 . The asymmetry S~ = J xs~dx = J x[s(x) — s(x)]dx was found to be —0.0003 ± 0.0011 
tending slightly negative, but consistent with zero. This differs from the published LO resul 
of -0.0027 ± 0.0013 which was obtained by a LO fit to combined NuTeV and CCFR data with 
a different set of PDFs. In this work we are only fitting the Nu TeV tables with CTEQ PDFs. 

With the LO baseline, an NLO cross sectioiPElfit using CTEQ6lJ21 nonstrange 

PDFs was then performed with the same parameterization as LO. S~ in this case was —0.0011 ± 
0.0014, again consistent with zero but tending negative at low x. The asymmetry is plotted on 
the right half of figure^ This method technically violates the Alterelli Parisi equations however, 
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Figure 1: Left: CTEQ5L LO fitted x(s-s) vs x Right: NLO CTEQ6M based fit x(s-s) vs x. 



by modifying pre-evolved PDFs by an arbitrary function of x. Another fit was performed with 
the same functional form for the strange and antistrange seas, this time defining s(x) and s(x) at 
a Qo of 1.3 GeV. The PDFs were then properly evolved in Q 2 with a version of the EVLCTEQ 
evolution code that allowed the strange and antistrange seas to be differen t The results 
obtained were very similar, with an asymmetry of —0.0013 ± 0.0013. 

Evolving the strange and antistrange seas properly brings the analysis closer to full QCD 
compliance, however the above scheme still violates the flavor sum rule by not requiring J s~ 
to be zero, and violates the momentum sum rule by not nor mali zing the strange sea relative to 
the nonstrange PDFs. Following the example of Olness et aff^ a fit was also performed using 
the following parameterization: 



~(x,Qo) = K + (l — x) a x 1 u(x,Qo) + d(x,Qo) 



s (x,Q ) 



s + (x) tanh 
s + + s~ 



K-(1-X) a X^~ 



xo 



(4) 
(5) 

(6) 



In this case, the flavor sum rule is enforced by choosing a crossing point xo such that / s~ 
is zero. Enforcing the momentum sum rule constraint was found not to affect the asymmetry. 
Results from fitting with this parameterization are shown in the left half of figure EJ The data 
prefers satisfying the flavor sum rule by forcing s~ to spike positive below an xo of 0.009, lower 
in x than any of the data table points. The asymmetry follows the same trend as the previous 
parameterization with an S~ of —0.0009 ± 0.0014. 

To examine whether a positive asymmetry is compa tible with the NuTeV dimuon data, an 
s~ shape similar to that found by the CTEQ collaboration^ was held fixed while s + was allowed 
to find the best minimum. The black curve shown in the right half of figureElis the result. This 
has a positive S~ of +.0016, less than one quarter of what is needed to resolve the NuTeV 
sin 2 6\v discrepancy. The best x 2 of this fit was high, 55/40 NDF. When s~ was then freed, the 
asymmetry quickly sought out the shape found on the left of figure |2] with a x 2 °f 37/37 NDF. 



4 Conclusion 



Several new fit results for the strange asymmetry using the NuTeV forward dimuon cross section 
table data have been presented. In general all the fits, LO or NLO, constraining sum rules or 
not, tended towards a negative strange asymmetry at low x that was still consistent with zero. 
Enforcing the flavor sum rule by constraining J s~{x)dx to be zero forced s~{x) to spike positive 




Figure 2: Left: xs (x) vs x. from fit with CTEQ s + ,s scheme. Right: Comparing asymmetry curves from fits 
with k — a parameterization (blue), CTEQ s + ,s~ parameterization, and a sample positive asymmetry 



below where there was data to constrain it. When the asymmetry was constrained to be positive, 
the x 2 increased significantly, making even a small positive asymmetry difficult to accomodate 
with our data. 

These results are still in the process of finalization, and it should be noted that the CTEQ 
collaboration has reached different conclusions from us in their preliminary global fits including 
the NuTeV tables. It is agreed however, that the NuTeV sin 2 6\y discrepancy is unlikely to be 
the product of an asymmetry in the strange and antistrange PDFs. NuTeV and CTEQ are 
continuing an active collaboration to try to resolve any remaining issues and to fully exploit the 
rich physics possibilities offered by the NuTeV dimuon data. 

References 

1. J. M. Conrad, M. H. Shaevitz and T. Bolton, Rev. Mod. Phys. 70, 1341 (1998) 

2. D. A. Harris et al. [NuTeV Collaboration], Nucl. Instrum. Meth. A 447, 377 (2000) 

3. S. J. Brodsky and B. Q. Ma, Phys. Lett. B 381, 317 (1996) 

4. S. Davidson, S. Forte, P. Gambino, N. Rius and A. Strumia, JHEP 0202, 037 (2002) 

5. G. P. Zeller et al. [NuTeV Collaboration], [Erratum-ibid. D 67, 119902 (2003)] 

6. S. Kretzer, F. Olness, J. Pumplin, D. Stump, W. K. Tung and M. H. Reno, 
|arXiv:hep-ph/03 12322} 

7. G. P. Zeller et al. [NuTeV Collaboration], Phys. Rev. Lett. 88, 091802 (2002) [Erratum- 
ibid. 90, 239902 (2003)] 

8. K. S. McFarland and S. O. Moch, arXiv:hep-ph/0306052| 

9. H. L. Lai et al. [CTEQ Collaboration], Eur. Phys. J. C 12, 375 (2000) 

10. J. Pumplin, D. R. Stump, J. Huston, H. L. Lai, P. Nadolsky and W. K. Tung, JHEP 0207, 
012 (2002) 

11. M. Gluck, S. Kretzer and E. Reya, Phys. Lett. B 380, 171 (1996) [Erratum-ibid. B 405, 
391 (1997)] 

12. S. Kretzer and I. Schienbein, Phys. Rev. D 58, 094035 (1998) 

13. M. Goncharov et al. [NuTeV Collaboration], Phys. Rev. D 64, 112006 (2001) 

14. S. Kretzer, D. Mason and F. Olness, Phys. Rev. D 65, 074010 (2002) 

15. B. Seligman PHD thesis, Columbia University, page 250 

16. T. Bolton, |ArXiv:hep-ex /9708014"l 

17. P. D. B. Collins and T. P. Spiller, J. Phys. G 11, 1289 (1985). 

18. W.K. Tung, private communication 

19. F. Olness et al., |arXiv:hep-ph/031232"3"l 



